All-proportional solid solution versus two-phase coexistence in the Ti–V alloy by first-principles phase field and SQS methods

The microstructures of the Ti–V alloy are studied by purely first-principles calculations without relying on any empirical or experimental parameter. The special quasirandom structure model is employed to treat the all-proportional solid solution \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β phase, while the first-principles phase field method or its variant is employed to treat the coexistence phases. The linearity of the calculated local free energy against the integer Ti\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_n$$\end{document}nV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_m$$\end{document}m composition in the cluster expansion method manifests a clear evidence of the solid solution behavior. From a detailed energy comparison, our results are consistent with the experimental fact that the Ti–V alloy is an all-proportional solid solution of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β phase at high temperatures and exhibits an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha +\beta$$\end{document}α+β coexistence at low temperatures. Moreover, it is found that mosaic-type microstructures may appear as a metastable phase, as observed by many experiments. The first-principles criterion for the all-proportional solid solution behavior presented in this paper is very general and can be applied to any other binary or multi-component alloys.


Methods
The FPPF method, the CPF model, and the SQS model are used in this study together with the first-principles method. A detail explanation of each method is given in what follows. The free energies obtained by these methods are compared for various compositions of the Ti-V alloy both at 0 K and 1300 K. The workflow adopting these methods is illustrated in Fig. 1 for easily understanding the present strategy. The purpose of this study is to compare the resulting free energies of three different FPPF, CPF and SQS models. In the FPPF and CPF simulations, we assume mixture of the hcp and bcc phases, while in the SQS simulation, we assume the bcc phase. Different models have different assumptions. Only after the final free energy comparison, we can judge which structure should realize in nature.
The FPPF method. The space is divided into fine grids, each of which forms a unit cell. We use a tetrahedron approximation in a bcc model (see the upper left panel of Fig. 1) in the cluster expansion method including atomic vacancies 24 . The unit cell is a √ 2 × √ 2 × 1 supercell of the conventional cubic unit cell and has four lattice points, which can be occupied by either Ti, V, and vacancy. Throughout this paper, we use the bcc lattice constant a = 4.425 Å obtained by the fully relaxed Ti 64 V 64 random configuration generated by the SQS model 28,29 . It is well known that the hcp α phase is energetically much lower than the beta phase when the V concentration is low. Therefore, for pure Ti, we use the result of our previous analysis for Ti64 alloy 21 . The concentration of the X element, φ X , represents the number of atoms inside the unit cell. That is, we assume that φ Ti ( φ V ) = 0-1, 1-2, 2-3, 3-4, and 4-5 correspond, respectively, to n (m) = 0, 1, 2, 3 , and 4 of Ti n V m (ideal configurations satisfy n + m = 4 , but we consider more general configurations 0 < n + m ≤ 6 as well). For the primitive cell calculations based on the tetrahedron approximation, the computational conditions are the same as those in our previous calculation of the Ti64 alloy 21 . (There, we have also compared the local free energy of the (low temperature) α (hcp) and (high temperature) β (bcc) phases of pure Ti, and confirmed that there is a phase transformation between these phases at certain temperature 21 .) In the FPPF simulation, it is not assumed that the structure is www.nature.com/scientificreports/ bcc. The alloy is a mixture of the hcp and bcc structures. Only in the SQS model, we assume the bcc structure, although the atomic positions inside the 4 × 4 × 4 supercell are fully relaxed. In order to determine the local free energy, we use potential renormalization theory 20,25 . Employing a supercell, we evaluate the total energy E j (r i ) where a central atom is displaced around the lattice point R and all other atoms are fixed at the surrounding lattice points. Taking account of the symmetry of the lattice, we displace the central atom inside the irreducible wedge of the Wigner-Seitz cell, and calculate the total energy deviation �E ij = E j (r i ) − E j (R) from the perfect crystal without displacement. Here, the suffix j denotes the species of the center atom and the suffix i denotes the i-th grid point around the lattice point R . Calculating E j (r i ) for many r i for each center species j = Ti or V, we evaluate the local partition function to determine the renormalized potential that depends on the temperature, which gives the contribution to the local free energy where C j is the concentration ratio of the center species j = Ti or V for a given configuration (Ti, V, TiV, TiV 3 , or Ti 3 V). Finally, by dividing the result by 4 20 , we obtain the desired renormalized correction ∆F in the basic tetrahedron. We add this F to the total energy E j (R) of the perfect crystal, in which all atoms were fixed at the lattice points, to obtain the local free energy at 1300 K. In general, explicit phonon calculations as in the Debye model are quite difficult due to the appearance of imaginary frequencies for unstable structures such as a pure BCC Ti or a B2 TiV at equiatomic composition, which may be less stable than a random phase. Therefore, potential renormalization theory has a distinct merit to treat the effect of off-lattice movements of atoms. The vibration entropy is included in the renormalization procedure under the framework of Einstein model, where the atoms can independently vibrate around each lattice point. It is very well known that the specific heat calculated by using the Einstein model is valid at least in high temperature regions including room temperature. The local free energy at 1300 K is calculated by taking the trace of the local partition function and eliminating the degrees of freedom at the length scale smaller than the lattice constant as described above. We use 2 × 2 × 2 bcc cubic supercell including 16 atoms to obtain the local free energy correction F for each Ti 4−n V n with integer n; see also 22 . As well as the total energy, the resulting local free energy becomes a stepwise function of the concentration, i.e., the number of atoms of each element inside one tetrahedron. The renormalization procedure is the same as our previous papers on the FPPF method [20][21][22][23] , which can automatically be done by using the Pipeline Pilot protocol 21 .
Then, we write the concentration, the chemical potential, and the mobility of element X as φ X , µ X , and M X , respectively. The flux is given by J X = −M X ∇µ X , and, from the continuity equation ∂φ X /∂t = −∇ · J X , we derive the Cahn-Hilliard equation where η represents the random fluctuation. If one regards the local free energy F as a continuous function of the concentration, and assumes that the chemical potential is given by µ X = ∂F/∂φ X , one obtains the standard Cahn-Hilliard equation 30,31 However, since the local free energy is a step function of the concentration, we have to take difference instead of differentiation with respect to the concentration to obtain the chemical potential. (Note that difference is a more natural definition for the chemical potential.) So, we calculate the chemical potential by 20 where ε X represents the gradient energy coefficient. The mobility and the gradient energy coefficient are assumed to be independent on the element. (At least for binary alloys, the interfacial energy is usually the same for the Ti and V sides.) They are absorbed by the time scale and the length scale, respectively, by rescaling the time and space coordinates. Therefore, they are dimensionless. The actual size can be determined by comparing the resulting patterns with experimental microstructures. Our model does not contain interfacial energy anisotropy or elastic energy contribution. However, the FPPF method is different from the standard PF model using the continuous local free energy. It is a distinct characteristic of the FPPF method that the faceted shape appears as a result without introducing any interfacial energy anisotropy or elastic energy contribution.
For the FPPF simulation, we use a 3D model composed of 40 × 40 × 40 meshes with the periodic boundary condition. The mesh space (dx) and the time step (dt) are set at 0.27 and 0.000125, respectively. The mobility M and the gradient energy coefficient ε , which are regarded as scaling parameters and can be set arbitrary, are set at 1 and 0.01, respectively. The amplitude of the random force is chosen as |η| = 0.5 . We assume initial fluctuations of the composition of ±c with c = 0.3 at 22 spot points with radius of 4 meshes. The FPPF simulation code is available on the website 32 , and input files for the FPPF simulation are written in Supplementary Information (SI).
The CPF model. In this study, the CPF model is specifically designed for the first time to simulate the solid solution phase of a binary alloy. We will see that the local (free) energy has almost a linear dependence on the composition. Therefore, the local (free) energy can be approximately replaced by a linear interpolation between the pure bcc Ti energy A and the pure bcc V energy B as www.nature.com/scientificreports/ where n stands for the V composition ( 0 ≤ n ≤ 4 ) of the Ti 4−n V n alloy (in a unit cell including one tetrahedron). In order to treat the coexistence of the α hcp Ti phase and the β bcc random Ti-V phase, we assume that a small fraction x ( 0 ≤ x ≤ 4 − n ) is the pure α hcp Ti phase. Then, the local (free) energy of the coexistence phase becomes at a given total Ti composition n. Here, C is the (free) energy of the α hcp Ti. (The local energies of A, B, and C at 0 K are changed to the local free energies by adding the renormalization correction F ren at 1300 K.) To obtain the total (free) energy, we have to add the interfacial energy for the interface between the α hcp Ti region and the β bcc random Ti-V region, which may be estimated as follows: The volume of the α phase compared to the total simulation cell volume is x/4, and therefore its surface area must be proportional to x 2/3 . However, the α phase does not appear as a single island but would have a lamellar or many-islets structure. In such cases, the number of islets would be proportional to x. This is of course an assumption, but anyway the power of x must be greater than 1 in the interfacial energy. Otherwise there is no minimum in the total (free) energy. Therefore, the total surface area of the α phase would be proportional to x 2/3 x = x 5/3 . Thus, we have to add Dx 5/3 to Eq. (6), where D is a proportional constant for the interfacial energy. This form of the interfacial energy is not essential at all and can be different. If we minimize this total (free) energy with respect to x at given (fixed) n, we have (5D This is the most suitable fraction of the pure α hcp phase in the coexistence. If the interfacial energy coefficient D is large, x m becomes small. On the other hand, if the interfacial energy coefficient D is small, x m becomes large and may be larger than 4 − n . In this case, since x cannot exceed 4 − n , x must be set equal to 4 − n . That is, the maximum value of x is 4 − n anyway. Here, x m = 4 − n only gives the boarder for the vanadium concentration n, where the pure α-Ti fraction x is treated as x = 4 − n or x = x m . We assume D = 0.05 eV corresponding to x m = 3.38 at 0 K (3.41 at 1300 K) as a typical value. Introducing the new variable φ = n + 0.5 , we define the continuous phase-field (CPF) model, which obeys the Cahn-Hilliard equation where F(n) is given by Eq. (6) with x = 4 − n for n > 4 − x m or with x = x m for n < 4 − x m . Note that F(n) does not explicitly include the interfacial energy term Dx 5/3 discussed above. The interfacial energy appears as the gradient energy proportional to ε in Eq. (8). Here, all the detailed settings of the CPF simulation (therefore the input files) are the same as for the FPPF simulation discussed above.
For the CPF simulation also, we use a 3D model composed of 40 × 40 × 40 meshes as well. All detailed settings of the simulation conditions are the same as those for the FPPF method mentioned above. The CPF simulation code is available on the website 32 , and its usage is written in Supplementary Information (SI).
First-principles method and the SQS model. For the first principles calculations, we use the planewave pseudopotential code CASTEP 33 . Using the Perdew-Burke-Ernzerhof (PBE) functional in the generalized gradient approximation (GGA) of DFT, we set the cut-off energy for plane waves at 500 eV in the on-the-fly generated ultrasoft pseudopotential. In order to simulate purely random configurations for various compositions, we employ the SQS model 19,28,29 . For this purpose, we assume 4 × 4 × 4 supercell of the conventional cubic unit cell, which includes 128 atoms. Then, the atomic geometries are fully relaxed with the fixed supercell frame using CASTEP. We have generated several samples and compared their total energies, but the total energy dependence on the sample was of the order of 10 meV/atom. Since the calculation of the renormalization at 1300 K is quite difficult for the random configurations in the 4 × 4 × 4 supercell, we approximately use the renormalization correction F ren interpolated according to the concentration. We also add the mixing entropy term �F mix = −TS mix = k B T{(4 − n) log[(4 − n)/4] + n log(n/4)} per four atoms. Then we can estimate the free energy at 1300 K.

Results and discussion
The FPPF method. The calculated energy values for the 5 integer compositions of Ti n V m with n + m = 4 in the tetrahedron approximation are listed in Table 1 together with the values for the α hcp phase 21 . The values inside the parentheses mean the local free energy including the effect of the potential renormalization at 1300 K. At most two interstitial atoms are considered inside the unit cell as well as the atomic vacancies. Its 2D map is shown in Fig. 2. The values for all 27 integer compositions of Ti n V m with 0 < n + m ≤ 6 are given in Table S1 of SI.
The most significant characteristic of the TiV alloy is the linear dependence of the local (free) energy on the V concentration, if we look at the integer Ti n V m ( n + m = 4 ) compositions in the tetrahedron 4-atom unit cell; see Fig 3a. This behavior has not been observed in any alloys in our previous studies; NiAl alloy 20 (see Fig. 3b),  23 . This linear dependence is a distinct characteristic of the low multiplephase-transformation barrier height or the seamless continuous phase transformations, and is a strong sign of the all-proportional solid solution. In the all-proportional solid solution, the two atoms can be mixed at any composition, and there is no stoichiometric composition. That is, there is no special meaning of the integer compositions like Ti 4 , Ti 3 V, Ti 2 V 2 , TiV 3 , and V 4 in the tetrahedron model. They are only some passing points in the all-proportional solid solution. This means that at any composition, the resulting local (free) energy should roughly coincide with the arithmetic weighted average between the two extremum compositions, Ti 4 and V 4 , which is indicated by a blue straight line in Fig. 3a. The real plots in between the two extremum compositions are slightly lower than the straight line, reflecting the additional energy gain by the mixing of the two atoms. So, roughly speaking, we can use this linear dependence of the local (free) energy of the complete A n B m ( n + m = 4 ) compositions as a criterion to select the candidates for the all-proportional solid solution. This criterion is rather qualitative and not very deterministic, but we can use it for the screening purpose. Three FPPF simulations are performed for 10 5 time steps where the microstructure do not change anymore. FPPF1 is the simulation based on the original tetrahedron model with four atoms in the unit cell. The local free energy values, which we used in the tetrahedron model for the FPPF1 simulation, are those of hcp Ti 4 , bcc     Figure 4p shows the initial V concentration in the case of 62.5 at% V. There are striking V precipitates in Fig. 4a of 6.25 at% V or Fig. 4i of 56.25 at% V and Ti precipitates in Fig. 4o of 93.75 at% V or Fig. 4k of 68.75 at% V, which are obviously an artifact of the tetrahedron approximation in a treatment of the Ti-V alloy. Therefore, we doubled the unit cell size to include 8 atoms and performed FPPF2. For the local free energy, we made a simple interpolation between the adjacent concentrations. The result is depicted by Fig. 5. Figure 5p is the initial V concentration in the case of 87.5 at% V. Now, all strange precipitates disappear, and mosaic and homogeneous patterns alternately appear at 6.25 + 12.5n and 12.5n at% V for n = 0, 1, 2, . . . , 7 . We further doubled the unit cell to include 16 atoms and performed FPPF3, whose result is depicted by Fig. 6. In this simulation, the patterns are homogeneous for all 6.25n at% V for n = 0, 1, 2, . . . , 16 . However, even in this case, a mosaic pattern appears at 90.625 at% V as shown in Fig. 6p. Although only one example is shown in this figure, all the half compositions of 3.125 + 6.25n at% V have a similar mosaic pattern. These mosaic-type microstructures www.nature.com/scientificreports/ have a possibility to appear at least in metastable phases. Indeed they very much resemble the experimentally observed twin-like microstructures of metastable phases 8-10 (compare, for example, Fig. 5 of Ref. 10 for 10 at% V with Fig. 4b for 12.5 at% V), although we cannot distinguish these phases with the metastable ω phase [34][35][36] within the present simple analysis. In the FPPF method, the final free energy is evaluated as the spatial integration of the local free energy distribution plus the interfacial energy. The calculated FPPF free energies including the interfacial energy F int = (ε/4) X=Ti,V [∇φ X (r)] 2 dr (additional factor 1/2 was multiplied to avoid double counting of Ti and V) at 0 K and at 1300 K are listed in Table 2. The values of the interfacial energy F int per four atoms are the order of 1 × 10 −4 eV for homogeneous patterns, about 0.08 eV for mosaic patterns, and 0.004-0.07 eV for patterns with faceted precipitates. The (free) energy obtained by FPPF3 are also written in the seventh and ninth columns of Table 3 for later comparison.
The SQS model. For a reliable judgement whether the target alloy is really the all-proportional solid solution, one should do more precise energy comparison. For this purpose, we performed 4 × 4 × 4 supercell calculations of the purely random configurations generated by using the SQS model 28,29 for variety of compositions. Then, the atomic geometries were relaxed with the fixed supercell frame. The resulting energy at 0 K for each www.nature.com/scientificreports/ composition is listed in the second column (0 K) of Table 3. This energy again has almost a linear dependence on the composition as seen in Fig. 3c. To obtain the free energy at 1300 K, we have to add the renormalization correction F ren and the mixing entropy term F mix . Then we obtain the free energy at 1300 K. All these values are listed in the third to fifth columns of Table 3.
The CPF model. The local (free) energy of the CPF model, Eq. (6), is shown in Fig. 7. The CPF simulations are performed for 10 5 time steps where the microstructure do not change anymore. The resulting patterns at 0 K on the z = 1 cross section are shown in Fig. 8. The initial V concentration is shown in Fig. 8p, which is the same as Fig. 4p. Although there is no clear microstructure, light shading patterns appear for all compositions. 1D profile along the blue allow of (o) 93.75 at% V is depicted in Fig. 9. The V-depletion region around y = 8 corresponds to the α phase and the other region corresponds to the β phase. Such a fluctuation in spacial V distributions can be considered as a characteristic of the α + β coexisting phase. Since we did not introduce any orientational order parameter 37 in our model for simplicity and clarity, there is no polycrystal-like domain structure as seen in the experimental patterns 5 . Besides this slight difference, the resulting almost homogenous pattern well resembles the experimental patterns; see, for example, Fig. 1(d)-(e) of Ref. 5 . There is no evidence www.nature.com/scientificreports/ of the ω phase 34-36 since it appears in metastable states only. There is also no twin-like structure as seen in the experimental patterns under quench, rolled or high pressure [8][9][10] . In the result of this CPF simulation, the interfacial energy calculated by (ε/2) [∇φ(r)] 2 dr is quite small (at most 1 × 10 −5 eV per four atoms) and negligible. The resulting total energy at 0 K and total free energy at 1300 K are summarized in the sixth and eighth columns of Table 3.
Energy comparison. Now we compare the resulting free energies. From Table 3, we immediately find that the FPPF free energy is higher than the CPF free energy and the free energy of the random configurations for  Table 3. Total free energy per four atoms of each composition in units of eV. For the purely random configurations, the energy at 0 K, the potential renormalization correction F ren , the mixing entropy F mix = −TS mix , and the free energy at 1300 K are listed in the second to fifth columns. For the results of the CPF and FPPF3 simulations, the energy at 0 K and the free energy at 1300 K are listed in the last four columns.  Figure 7. Continuous local energy used for the CPF model at 0 K. Continuous local free energy at 1300 K looks the same, because the potential renormalization correction F ren is very small. www.nature.com/scientificreports/ all V compositions at both 0 K and 1300 K, although the difference is not very large but only 0.02-0.59 eV per four atoms. Therefore, the mosaic-type microstructures shown in Fig. 6p may be realized at least as the metastable phases as observed in many experiments [8][9][10] . Next we compare the CPF free energy with the free energy of the random configurations. At 0 K, the former is lower than the latter for most V concentrations, although the former is the same as or slightly higher than the latter for 62.5-87.5 at% V by 0.01-0.02 eV per four atoms. However, the CPF patterns ( Fig. 8) are almost homogeneous and cannot be distinguished from the purely random configurations. Therefore, we conclude that the almost homogeneous α + β coexistence phase appears at 0 K. On the other hand, at 1300 K, the CPF free energy is higher than the free energy of the random configurations for most V concentrations, although the former is lower than the latter for 6.25-31.25 at% V by 0.05-0.33 eV per four atoms. It is true that there is an error in the free energy of the random configurations due to the approximation introduced in the renormalization calculation. However, again, the CPF patterns are almost homogeneous and cannot be distinguished from the purely random configurations. Therefore, we can conclude that the almost homogeneous β random phase appears at 1300 K. Overall, our results are consistent with the experimental fact that the Ti-V alloy forms the α + β two-phase coexistence at low temperatures, and the all-proportional solid solution at high temperatures.

Concluding remarks
In this paper, we have studied the phase behavior of the Ti-V alloy from purely first-principles without relying on any empirical or experimental parameter. We used the first-principles phase field (FPPF) method and its continuous variant (CPF model) to treat the coexisting phases, and the SQS model to treat the all-proportional solid solution β phase. Within the accuracy of the approximation used in the present analysis, the resulting energy comparison successfully derives the experimental fact at phase equilibrium that the Ti-V alloy is an all-proportional solid solution of the β phase at high temperatures and exhibits an α + β coexistence at low temperatures. Moreover, our result predicts that the mosaic-type microstructures shown in Fig. 6p may be realized at least as the metastable phases as observed by many experiments [8][9][10] . We also found that the linear energy dependence of the local free energy on the integer A n B m ( n + m = 4 ) compositions in the cluster expansion method can be used as a criterion whether the alloy exhibit a solid solution behavior. Although this criterion is rather qualitative and not very deterministic, it can be used for the screening purpose. We summarize the present general strategy to identify the all-proportional solid solution as follows. First, perform the first-principles total energy calculation of the primitive cell in the first-principles phase field (FPPF) method. Second, perform a very simple screening by checking the linearity of the energy calculated in this method. If this screening is passed, perform first-principles supercell calculations for purely random configurations of the candidate alloy using the SQS model. Then, construct the CPF model if necessary and perform the FPPF and/or CPF simulation. Finally, comparing the resulting energies (or free energies for finite temperatures), identify the most stable microstructure to be either a two-phase coexistence obtained by the FPPF method or a purely random configuration everywhere.
To demonstrate the validity of this strategy, let us briefly examine the case of the NiAl alloy 20 , which exhibits clear two-phase ( γ /γ ′ ) coexistence. First of all, the local free energy of the stoichiometric compositions Ni n Al m ( n + m = 4 ) (see Fig. 2 of Ref. 20 ) does not change linearly in n; the deviation from the linear interpolation is quite large (2.8 eV at Ni 2 Al 2 ); see Fig. 3c. Therefore, from the first screening, this alloy is clearly excluded. Furthermore, our 3 × 3 × 3 cubic fcc supercell calculation of the SQS-based purely random configurations gives the total energies of −23.54 eV and −26.58 eV, respectively for 50 and 83.3 at% Ni, and the corresponding values for the completely ordered B2 phase, and the L1 2 ( −26.26 eV) and pure Ni fcc ( −27.72 eV) coexistence phase are −24.14 eV and −26.70 eV (we added the coexistence interfacial energy of 0.05 eV estimated by the original FPPF method to the arithmetic weighted average of the L1 2 and pure Ni fcc phases), respectively, which indicates that the B2 phase at 50 at% Ni and the two-phase ( γ /γ ′ ) coexistence at 83.3 at% Ni are both energetically stable at 0 K. This situation does not change even if we add the mixing entropy of − 0.13 eV and − 0.09 eV, respectively, to the random configurations at 1300 K. Then, assuming that the effect of the potential renormalization at 1300 K is the same both for the uniform random phase and for the coexistence phase, we are led to the conclusion that the coexistence phase is more stable at 1300 K also. This strategy is very general and applicable not only to binary alloys but also to ternary and multi-element alloys. For the FPPF method, we have already succeeded in applying it to ternary alloys like Ti64 21 . To extend the CPF model to ternary alloys is www.nature.com/scientificreports/ left for future study, but the strategy is straightforward although it depends on case by case. Application of SQS model to ternary alloys is also possible. So, although the application of the present method to ternary alloys is left for future study, there is no essential difficulty in doing it. We believe that the present method is useful for determining the phase morphology of alloys.

Data availability
All data generated and analysed during this study are included in the published article and its Supplementary Information.